Critical phase in non-conserving zero-range processes and equilibrium networks 
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Zero-range processes, in which particles hop between sites on a lattice, are closely related to 
equilibrium networks, in which rewiring of links take place. Both systems exhibit a condensation 
transition for appropriate choices of the dynamical rules. The transition results in a macroscopically 
occupied site for zero-range processes and a macroscopically connected node for networks. Criticality, 
characterized by a scale-free distribution, is obtained only at the transition point. This is in contrast 
with the widespread scale-free real-life networks. Here we propose a generalization of these models 
whereby criticality is obtained throughout an entire phase, and the scale-free distribution does not 
depend on any fine-tuned parameter. 



PACS numbers: 89.75.-k, 05.70.Ln, 05.40.-a 

Many driven, non-equilibrium models, reach a critical 
or scale invariant steady state only when their dynami- 
cal parameters are fine-tuned to reach a phase transition 
point. Examples include a wide range of systems such 
as jamming in traffic [J, coalescence in granular gases 
0, gelation in networks and wealth condensation in 
macroeconomics Q. In all these systems one has a con- 
densation phase transition which we shall discuss in de- 
tail below. In other non-equilibrium models, for example 
in driven lattice gases and sandpile models 0, , it has 
been argued that scale invariance and power-law distri- 
butions are generic, or at least one may have scale-free 
distributions across wide regions of the parameter space 
rather than just at critical points. This phenomenon has 
been termed self organized criticality. 

In recent years considerable attention has been given 
to the study of real-life networks. Networks, defined as 
collections of nodes connected by links, are found in many 
fields of study, ranging from molecular biology to social 
communities and the Internet With each node one 

associates a degree k which is the number of links con- 
nected to it. In general links may be directed or they may 
carry a weight, however for our purposes we do not con- 
sider such features. It has been observed that very often 
real-life networks are characterized by a degree distribu- 
tion p{k) which decays algebraically for large fc 0, 
These networks, termed scale-free networks, are indeed 
critical, suggesting the existence of a mechanism which 
drives them to this state. Subsequently, dynamical pro- 
cesses for growing networks have been proposed which 
result in a critical distribution for a wide range of the dy- 
namical parameters. In these processes nodes and links 
are continually added to the network with some prede- 
termined rates 0, Isj. [To[ . On the other hand "equi- 
librium" networks [3, 1131 j whose dynamics constitutes 
rewiring processes with a fixed number of nodes, exhibit 
a critical distribution only at a critical point. This tran- 
sition corresponds to condensation (also referred to as 
gelation) where a single node captures a finite fraction 
of the links. Note that although termed equilibrium net- 
works, their dynamics does not always obey detailed bal- 



ance, thus their steady state may not always be a thermal 
equilibrium one. 

Instructive insight has been gained into the condensa- 
tion transition throu gh t he analysis of simple interacting 
particle systems 0,^3 0|. These systems form fun- 
damental models which may be mapped onto particular 
applications. For example, the zero-range process (ZRP) 
fl6i | is a particularly simple and exactly soluble model 
in which each site of a lattice may contain an integer 
number of particles and particles hop to a neighbor- 
ing site with rate u{n^j). This model is closely related to 
the equilibrium networks discussed above which undergo 
rewiring dynamics "T^l; in the following we will exploit 
this relationship. 

Condensation in the ZRP will occur, for example, when 
u{n) decays to some finite, large-n asymptotic value (3, 
as u{n) ^ P[l + b/n) with b > 2. The transition is simply 
understood by considering p(n) the steady-state proba- 
bility that a site contains n particles. For a subcritical 
density of particles one finds that p{n) decays exponen- 
tially with decay length dependent on the conserved par- 
ticle density p = N / L where N is the number of particles 
and L is the number of sites. This p{n) describes the low 
density, fluid phase. As the density is increased the de- 
cay length increases until at the critical density, pc, it di- 
verges and one has a power-law distribution p{n) ^ , 
thus a critical fluid. Above pc, in addition to the power 
law, a piece of p(n) emerges centered about n = L(p—pc); 
this piece represents the condensate [l^. Thus the su- 
percritical phase corresponds to a critical fluid coexisting 
with a condensate. Only at criticality does one have a 
pure power-law distribution. 

In the present work we investigate how a critical phase 
may emerge in processes such as ZRPs and network 
models. We shall see that the introduction of non- 
conservation in an appropriate fashion, modifies the con- 
densate phase into a scale-free phase by effectively sup- 
pressing the condensate and leaving the critical fluid. In 
this generalization of the ZRP particles are created at 
all sites at constant rate but are removed at a rate that 
depends non-linearly on the occupation number. Equiv- 



alently in the network context links are created and de- 
stroyed. We begin by elucidating the mechanism for the 
generation of a critical phase within a generalized ZRP; 
later we will define an equilibrium network where this 
mechanism is also manifested. 

Consider a lattice of L sites upon which reside a num- 
ber of particles. With rate u{n^) (probability per unit 
time) which depends on the occupation of site /x, a 
particle is transferred from site ji to another site. We 
consider a fully connected geometry, where the destina- 
tion site is chosen randomly from the other L — \ sites. 
In addition to the hopping dynamics, particles are added 
to site /X with a constant rate c, or removed with a rate 
a(n^), which increases with the site occupation n^. Thus, 
our choice for the dynamical rates is given by 

n, m ^ n — 1, m + 1 with rate u{n) = (l + ■^) 

. + 1 with rate c = {j^Y (1) 



n 
n 



— 1 with rate a{n) = 



where b, s and k are positive parameters, and 9{n) is the 
usual Heaviside step function. The dynamical rates are 
conveniently implemented by using a random-sequential 
updating scheme, whereby at each time step a site fi is 
chosen at random, and a hop, annihilation or creation 
event may occur with relative probabilities given by the 
rates With b < 2 the system is always found in a 
sub-critical phase, where the particle number distribu- 
tion is exponential. We therefore restrict our discussion 
hereafter to the case b > 2. 

In Fig. 1 we compare the particle number distribution 
p{n) of our model with that of the ZRP with conserving 
dynamics. It is clearly seen that for this choice of param- 
eters, the creation/annihilation dynamics can selectively 
destroy any condensate and sustain the power-law dis- 
tribution, corresponding to the critical fluid. In what 
follows we analyze the model showing that this feature 
holds for an entire region in the parameter space. 

The steady state of the model is fully described by the 
probability distribution P{ni, n2, . . . ,71^) over all possi- 
ble configurations. In contrast to the conserving ZRP 
[TQ] . the steady-state distribution of the model (Q) does 
not factorize generally. However, we make the mean field 
approximation that the steady state distribution does 

factorize, i.e. P(ni, n2, . . • , n^) Y[f=iP(.^i)- -'-'^^ 
the fully-connected geometry we expect this approxima- 
tion to become exact in the limit L — > oo. 

Using this approximation, the steady-state master 
equation is given by 

= [u{n + l) + a{n + l)]p{n+l)-[X + c]p{n) (2) 
- {[u{n) + a{n)]p{n) - [A + c]p{n - l)}0{n). 



Here the current A is given by 

oo 

A = u{n)p{n) 



(3) 




FIG. 1: Steady state distribution of the non-conserving ZRP 
(black solid line) and the conserving ZRP (grey dotted line). 
The data is from simulations run on a system of size L = 10*, 
with b = 2.6, fc = 3 and s = 1.96. In the conserving model 
the particle density was set to p = A > pc. The peak at 
high occupation number, which exists only in the conserving 
model, corresponds to the condensate. 



From (O it follows that 



p{n) 



(A + c)" 



p(0) 



(4) 



Note that this is not a closed solution as A depends on 
p{n). The values of A and p(0) should be set such that 
both the normalization condition 1 = X]n=oP(") ^-"^^ 
creation/annihilation balance condition 



a{n)p{n) 



(5) 



are obeyed. We now identify the three phases of the 
model by determining the asymptotic, large L behaviors 
of p{n) and A which satisfy (PJ and ©. The emergent 
phase diagram is summarised in Fig. 2. Deferring details 
to a later publication, we find the following results: 

Low-density phase, s > k — Rewriting Q as 
^fc-s _ ^ n^p(n) implies p{n) is a rapidly decreasing 
function of n, and the steady state density p is ^ 1. 
Thus, p{l) p ^ L^~^ and in the thermodynamic limit 
the density goes to zero. 

For s <k the sum in (jSJ is controlled by the behavior 
of p{n) at large n. We find the following regimes: 

High-density phase, s < k/{k + 1) — Here 



p{n) 



n ^ exp 



,k+l 



{k + \)L^ 



(6) 



thus p{n) is strongly peaked at n ^ i^"*/*-'. This is the 
high density phase, where all sites are highly occupied. 
Note that the mean number of particles in the system, 
N Lp—'^/^ is super-extensive. 

Critical phase, k > s > k/{k+l) — In this phase the 
system relaxes to the critical density, and p{n) takes an 
algebraic form. However, for large but finite systems two 




FIG. 2: Typical phase diagram of both ZRP and network 
models in the k-s plane with b fixed, k,s,b are defined in Q. 



sub-phases are observed, distinguished by the finite-size 
corrections to the dominant power law. 
For A: > s > kb/{k + 1) 



p{n) 



n ^ exp 



n 

1? 



(7) 



where x = {k ~ s) /{k ~ b + 1). This cuts off the power 
law at n '--^ L^. We refer to this as critical sub- phase (a). 
For kb/{k + 1) > s > k/{k + 1) 
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exp 
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^k+l 



{k + l)L^ 



(8) 



where d = b — s{k + l)/k. Here, on top of the algebraic 
part, p{n) is weakly peaked at n L''/(''+^'(lnL)^/(''+^). 
This peak will diminish as L — > oo. We refer to this as 
critical sub-phase (b). 

In Fig. 3 we present typical data obtained from nu- 
merical simulations in the different phases and compare 
with theoretical curves of p{n). We found that in all 
phases, starting from random initial configurations of 
various densities, the system relaxes towards its expected 
steady state value. However, for the low-density phase 
the time scales for full relaxation were prohibitive and we 
do not present steady state data for this phase. In Fig. 3 
we also provide data for a one-dimensional (Id) variant 
of the model, where sites are arranged in a Id array, and 
particles are allowed to hop only to the right neighbor of 
the departure site. For the Id geometry the mean-field 
approximation is not expected to be exact even in the 
limit L — > 00. Nevertheless, we find numerically that the 
three phases discussed above exist also in the Id model. 

We now apply the approach discussed above to equi- 
librium networks. To make the analogy with the ZRP, 
one identifies a site of the ZRP and its occupation num- 
ber with a node in the network and its degree, respec- 
tively. We define a network model which incorporates 
both rewiring and creation/annihilation dynamics, and 
show how a proper choice of rates leads to the existence 
of a critical phase, much like that of the non-conserving 




10000 



FIG. 3: Steady state distributions from simulations of the 
ZRP model on a fully-connected lattice and a Id lattice 
(O), compared with the theoretical asymptotic curves. Here 
L = 5000 and 6 = 2.6, fe = 3. Dashed curve is critical sub- 
phase (a) (s = 2); dotted curve is critical sub-phase (b) (s = 
1.2); full line is high-density phase (s = 0.4). 



ZRP, within which networks are scale free. Due to the 
introduction of annihilation of links, there is no simple 
mapping from the ZRP to the network model since this 
would require keeping track of pairs of linked particles 
in the ZRP. However the two systems are closely related 
and as we shall see share the same phase diagram. 

We consider a network of L nodes which are linked 
together by an integer number, iV/2, of undirected links 
{N is the number of particles in the corresponding ZRP). 
With rate u(n^) = I + b/n^ one of the links is dis- 
connected from node /i and is rewired to another ran- 
domly chosen node. This does not change the number of 
links in the network. With rate a(n^) = (n^/L)*^ one of 
the links connected to node is removed from the net- 
work. In addition, a new link is created between node 
fi and other randomly chosen node at a constant rate 
c = 1/L'^. Again the dynamics is conveniently imple- 
mented by choosing a node fi randomly at each time step, 
and changing the wiring with probabilities constructed 
from the relevant rates. The mean-field master equation 
for the network model differs slightly from that of the 
ZRP (0) , and is given by 



= [a(n + 1) +M(n+ 1) + A(n+ I)]p(n+ I) 

— [a{n) + u{n) + A(n)]p{n)9{n) 

- [X + 2c]p{n) + [X + 2c]p{ji~l)ein) . 

Here A = ^ u{n)p{n) as before, and 

oo 



(9) 



1=1 



(10) 



The steady state solution to ^ is given by 
(A + 2c)" 



p{n) 



rirLi [A("^) a(m) -|- u{m)] 



p(0) 



(11) 
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FIG. 4: Steady state probability distributions from simula- 
tions of the network model in the critical phase (sub-phase 
(a)). Here b = 2.6, fc = 3 and s = 2, with L = 1000 (circles), 
L = 2000 (squares) and L = 4000 (triangles). 



The main difference between this result and that of the 
ZRP (0J lies in A(n) and its dependence on the to- 
tal number of links in the system. This complicates 
slightly the analysis, however all phases persist includ- 
ing the critical phase. Thus the phase diagram in Fig. 2 
also describes the network model: the low-density and 
high-density phases characterize extremely sparse and ex- 
tremely dense networks; the critical fluid corresponds to 
a scale- free phase in the thermodynamic limit. Examples 
of the degree distribution in the critical phase are given in 
Fig. 4, where we present data obtained from simulations 
of systems of increasing sizes. The power-law regime in- 
creases with system size. 

Further interesting observations are that the weak 
peak of p{m) in critical sub-phase (b) may correspond 
to a number of highly connected nodes but is distinct 



from a condensate phase. This would correspond to a 
large number of hubs in the network. Also, we note that 
in the network model the suppression of events which 
link nodes to themselves has a considerable effect on the 
results , by introducing a new cut-off into the proba- 
bility distribution which can become the dominant scale. 
A full analysis will be published elsewhere. 

Our main interest in the systems studied lies in the 
emergence of a critical phase which we have shown ex- 
ists for annihilation and creation indices fc, s in the range 
k > s > k/{k + 1). We conclude by comparing the crit- 
ical phases we have identified to the critical points of 
the corresponding conserving ZRP and network models. 
In the latter, the average particle/link density p is an 
external parameter. A power-law distribution of the oc- 
cupation number/degree is only obtained at p — pc- In 
contrast, in the non-conserving models we have studied, 
the steady-state density is set by the dynamics to be pc 
and a power-law distribution is obtained throughout the 
critical phase. In models exhibiting self-organized criti- 
cality [l7| , a critical phase is typically obtained only when 
the driving rate of the system vanishes with the system 
size |T^ . in order to ensure relaxation between stimuli. 
In comparison, in the present work the creation rate c in 
e.g. QJ, which could be thought of as a driving rate for 
the system, vanishes in the large L limit whereas u the 
hopping/rewiring rate does not vanish. Thus, there is a 
separation of timescales in the dynamical processes. On 
the other hand, there are no avalanches or underlying ab- 
sorbing states which are features usually associated with 
self organized criticality. 
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